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Abstract 

A simple method is used to extrapolate the form factor for the process 7*7* — > 
7T° when two virtual photons produce a 7r° from the region of large spacelike 
photon virtualities to the experimentally accessible case when one of the photons 
is real. The method is also applied to the study of the axial anomaly and uj — 
p — it coupling and to the pion charge and charge radius. The results compare 
favourably with the experimental measurements. 



a: 

j§ : 1 Introduction 



The study of the process 7*7* — ► tt° and the related form factors when two virtual 
photons produce a neutral pion provides a simple, interesting and clean test of QCD. 
In contrast to the pion electromagnetic form factor it is not asymptotically supressed 
by being multplied by the strong coupling constant a s . In the lowest order in a s , 
perturbative QCD predicts that [[]]] 

F^AM) = A if ^ 2 dx (L1) 

11 3 Jo xQ 2 + xq 2 

where (j) n (x) is the phenomenological pion distribution amplitude and x, x = 1 — x 
are the fractions of the pion light cone momenta carried by the quarks in the region 
where both photon virtualities q 2 = —q\ and Q 2 = —q\ are large. 

The most experimentally favourable situation is when one of the photons is real, 
q 2 = and this has led to the measurement of the form factor F 11 * n o(q 2 = 0, Q 2 ) by 
the CELLO [§ and CLIO [§ groups. 

Recently, Radyushkin and Ruskov [Q] have undertaken to obtain the transition form 
factor F 77 * 7r o(g 2 = 0, Q 2 ) independently of the specific form of the pion distribution 
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amplitude. They have rather first evaluated the perturbative and non-perturbative con- 
tributions to the amplitude F~* 7 *_^o (q 2 , Q 2 ) in the deep Euclidean region and have the 
extrapolated to the mass-shell of one of the photons in a QCD sum rule [|5p approach. 
It is the purpose of this note to undertake this extrapolation in a more straightforward 
and much less model dependant way. 

The task of extrapolating QCD prodictions from the deep asymptotic region to low 
energies is not a trivial one, an effective way is the use of QCD sum rules 0] in which 
use is made of the analyticity properties of an amplitude to write dispersion integrals 
the low energy part of which is saturated by known resonances and where QCD is 
employed to evaluate the asymptotic part. It is in general hard to assess the quality of 
the approximations made as resonance saturation should not be expected to hold over 
too large an interval and deviations from QCD could become important far away from 
the Euclidean region. 

A simple model independant method of extrapolation is proposed in this note in 
which the integration kernel in the dispersion integral is modified by adding to it poly- 
nomials of low order N . The polynomials are chosen such as to annihilate the kernel 
as well as its derivatives up to order N in the region where both resonance saturation 
and QCD are least reliable. 

The method is applied to obtain the form factor F 77 . 7r o(g 2 = 0, Q 2 ). It is shown that 
extrapolation of the QCD expression of Radyushkin and Ruskov [0] with the values 
they use for the damping parameters, which they deduce from stability considerations, 
as well as the values commonly used for the QCD non-perturbative vacuum conden- 
satesQ lead to values of F 77 * 7r o(g 2 = 0, Q 2 ) in good agreement with the experimental 
data. 

Another application of the method is to the anomaly relation [^p F 777r o = , 
where /„. = 92.4MeV is the pion decay constant, which holds when both photons are 
on their mass-shell. In particular we evaluate corrections to the value of the to — p — n 
coupling constant given by vector meson dominance. 

The effectiveness and stability of the method is finally demonstrated by calculating 
the charge and charge radius of the pion from the asymptotic behaviour of the pion's 
electromagnetic form factor obtained some time ago by Ioffe and Smilga [|7p and by 
Nesterenko and Radyushkin []8||. Agreement with experiment is good in all cases. 

2 The 77*7T° form factor 

The electroproduction of a neutral pion in the reaction 7* — > 771-° is experimentally 
accessible [[p, [|3p. The corresponding form factor is defined when both photons are 
off their mass-shell. 

\p2ie^ a pq*q%F rr ^(q{,q 2 2 ) = 4vr dx exp(igix)(7r(p) | TV^{x)V u {0) | 0) (2.1) 



The result turn out to be insensitive to the exact values of these. 
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where is the electromagnetic current of light quarks divided by the unit charge 

2 

Radyushkin and Ruskov [Q] calculate the QCD expression of the form factor which 
is valid when both t = q\ and Q 2 = — q\ are in the deep euclidean region 



F^(t, Q 2 ) = yfaM&P^t, Q 2 ) = F p (t, Q 2 ) + c. 2 t + Co + j + ^ + 



with the perturbative part 



F p (t,Q 2 
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ds exp 
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(2.2) 



(2.3) 



M 2 is the Borel damping mass and R is the duality radius in the axial-vector 
channel. Both are determined by stability considerations []7|], [||], the values M 2 ~ 
.7GeV 2 are used. The coefficients 
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are expressed in terms of the QCD gluon and 4-quark condansates. There is no con- 
sensus on the value of the 4-quark condensate, in particular the vacuum saturation 
approximation ((qq) 2 ) = (qq) 2 is likely to be badly violated Fortunately the con- 
tribution of the non-perturbative condensates to our sum rules is quite small which 
renders our results insensitive to their exact values. 

As stated in the introduction, our purpose is to extrapolate expression ( [2.2D to t = 
and compare with the experiemtal measurment of F(0, Q 2 ) done by CELLO [Q] and 
CLIO [g]. 

At low values of t, F(t, Q 2 ) is dominated by the p— meson pole 

F (t,Q 2 ) = jf^\ + --- (2-5) 
(t - m 2 ) 

Furthermore, F(t, Q 2 ) is an analytic function in the complex t— plane with a cut on 
the positive real axis running from threshold t = Am 2 to infinity. Consider the integral 

1 



F(0,Q< 



2ni 



dtK N {t)F{t,Q z 



(2.6) 
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where c is the contour consisting of a circle of large radius R and two straight lines 
parallel to and lying just above and just below the cut. 



1 1 N / t \ n 

ra=0 v 7 

and the coefficients a n are up to now arbitrary constants. If R is large enough, the 



QCD expression ( PT2| ) is expected to hold on the circle except in the region close to the 
positive real axis. On the cut and for low values of t, p— pole dominance, expressed 
by eq. ( |2.5[ ) adequatly represents the amplitude. We now choose the coefficients a n 
defining the kernel eq. ( J2.7| ) such that K N {t) as well as its derivatives up to order 
N vanish at t = R in order to annihilate the integrand in the integration region not 
covered by vector meson dominance and where the QCD expression cannot be trusted. 
This is achieved by taking J2n=o a n ("1)™ t0 coincide with the Taylor expansion of | 



about t = R truncated at N. The integral fl2.q ) then yields 

CLqC2 CI1C4 



F(0, Q 2 ) = -g(Q 2 )K N (m 2 p ) + ^ j> dt K N (t)F p (t, Q 2 ) + c 



R R 2 

(2.8) 



The duality radius in the p— meson channel R is taken to be R ~ 1.5GeV 2 [§]. In 
order to determine F(0, Q 2 ) we take N — 1 and N = 2 successively in eq. Q2.8P and 



perform the integral over the circle of radius R in the complex t— plane. The result is 
shown in fig. 1 and agreement with the data is seen to be quite good. 

The deep inelastic limit Q 2 R is readily obtained, in this limit the pole of 
F p (t,Q 2 ) lies outside the integration contour and only the l/t term in K N {t) con- 
tributes so that 

The fact that the r.h.s. above should not depend on iV implies that Q 2 g(Q 2 ) 
as expected and that 

Q 2 F(0,Q 2 ) M 2 (l-exp (-f^)) 



8vr 2 / 2 Q 2 -oo 8tt 2 / 2 



Q 



.75 (2.10) 



The above result can be compared to values obtained from eq. ( |1.1[ ) if one uses 
asymptotic forms of the pion distribution amplitude as given by Brodsky and Lepage 
[gg] and by Chernyak and Zhitnitsky tf\A§ 
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Q 2 F BL (0,Q 2 ) 

8n 2 f 2 
Q 2 F cz (0,Q 2 ) 5 

8n 2 f 2 qCoo 3 



(2.11) 
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The limiting value given by eq. ( |2.10[ ), together with the bound at Q 2 = imposed 
by the anomaly relation F~*~* 7r o(0, 0) = J can be accomodated by the interpola- 
tion formula 



Q 2 F rr7T o(Q 2 ) _ Q 2 1 



ln 2 f 2 Q 2 + p 2 8rr 2 f 2 



(2.12) 



with fx 2 = M 2 (l — exp (-jfa)) = .5GeV 2 

Expression ( |2.12| ) yields values of the form factor very close to the ones obtained 
analytically from eq. |2~8| . 



3 The Anomaly Relation and u — p — n coupling 



The anomaly relation [|5p determines the form factor when both photons are on their 
mass-shell 



F(0,0) = v^tt/^XO, 0) 



(3.1) 



For the symmetric case q\ = q\ = t, the asymptotic expression of Radyushkin and 
Ruskov [ffll takes the form 



F QCD (t) 
with 
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As before we start from the integral 



dtK N (t)F(t) 



with Kjq{t) defined by eq. [2~77| . In this case, the amplitude F(t) possesses a double as 
well as a single pole at t = m 2 because the p and u mesons are degenerate in mass, 
i.e. in the low energy region 



F(t) 



f 



+ 



+ 



(3.4) 
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One obtains then 



F(0) = -fK' N (m*) - gK N (m 2 p ) + ± j> dtK N (t)F p (t) 



a c 2 aic 4 a 2 c 6 
(3.5) 



with the perturbative part F p (t) given by the first term in the r.h.s. of eq. Q3.2P 



The residue g which describes the transition 7r+vector meson^continuum is un- 
known. The residue / is related to G^^, the cuprc coupling constant, 

f = Sir 2 f 7T ——^-G UJ p 7T (3.6) 
9 P 9u 

where denote the coupling of the vector currents to the corresponding vector mesons. 
Taking N = and iV = 1 in eq. [33| with F(0) = 1 eliminates g and yields 



G up)t = 12.2G6V- 1 (3.7) 

a value close to the one given by vector meson dominance G wp7r = HGeV -1 . A check 
of the stability of the method is provided by the taking N = 2 with the values obtained 
for / and g inserted in the r.h.s. of eq. |33]to evaluate F(0), this yields^ 

F(0) = .96 (3.8) 



4 Pion charge and pion charge radius 

In this paragraph the effectiveness and stability of the method used above is tested in 
the calculation of the pion charge and charge radius. The electromagnetic form factor 
of the pion is 

(^ + (P 2 ) I V™ | n+( Pl )) = F(t)( Pl+P2 )p (4.1) 

and 



F(0) = 1 (4.2) 

is the charge of the pion. 

The asymptotic behaviour of F(t) was obtained some time ago by Ioffe and Smilga 
[0] and by Nesterenko and Radyushkin pSp 

F «c» (i) = _ mff _^ + m + SL + *_ i + l_SL t + ... (43) 



2 This corresponds to the factorisation value ((qq) 2 ) = (qq) 2 - The case ((qq) 2 ) = 3(qq) 2 , e.g., 
yields F(0) = .98. 



6 



the first term is asymptotically leading [JTp and 



F P (t) 



Ro exp (— jgs) £ 2 <Ar exp (— jfs) £ 3 cfe 



o 



(t - 2a;) 2 y (t - 2x) 3 



(4.4) 



1 t a s G\ 104vr 2 

c *=24^<— >' c *=8i^«w)> 

i?o — 1.2GeV 2 , M 2 ~ lGeV 2 have been obtained from stability considerations. 
As before consider 

F(0) = ^ J dtK N (t)F{t) 

K N (t) = J-^£a n (£) or 

F(0) = f p G p ^K N {m^ + ^-jdtK N {t)F QCD {t) (4.5) 
Using standard values for the condensates yields 



(4.6) 



F(0) = 


.95 


N-- 


= 


F(0) = 


1.15 


N 


= 1 


F(0) = 


1.15 


N 


= 2 



in reasonably good agreement with experiment. 
The derivative F'(0) yields the charge radius 



F'^) = ^-JdtK N {t)F{t) (4.7) 



where now 

A? 



I 1 / t \ n 

k ^) = W-W^{r) (4 ' 8) 

and once more the a n 's are chosen so as to annihilate and it s derivatives up to 
order N at t = R. 

F'(0) = f p G p7r7T K N (m 2 p ) + dtK N (t)F QCD (t) (4.9) 



the result is 



F'(0) = 1.7GeV 2 N = 
F'(0) = 1.93GeV 2 N = 1 

F'(0) = 2.05GeV 2 N = 2 (1) 
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corresponding to a charge radius 



( r 2) = (.47 ± .04) /m 2 (4.11) 

which compares quite well with the experimental value (r 2 ) = AAfm 2 . 

In the examples treated above the contribution of the condensates is quite small 
and for the chosen values of iV unknown higher condensates do not contribute. 

I thank the Academie de Versailles for its generous support and the IPN, Orsay for 
their kind hospitality. 
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